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ABSTRACT 


This thesis examines two problems associated with 
crystalline defects in the materials Zr and Zr3Al, by means 
of a molecular dynamics computer simulation. 

In the first problem, the relaxed configuration around a 
Single vacancy, a di-vacancy and a tri-vacancy 1S computed by 
introducing the defect into a perfect lattice and dynamically 
relaxing the crystal until it attains equilibrium. The 
results are in qualitative agreement with "static" studies 
but show clearly that first and second neighbors have 
different relaxed states depending on their configuration in 
enhewrmitial crystal. lattice. 

The recombination studies show that a vacancy- 
interstitial pair quickly relaxes to the perfect crystal if 


the interstitial is within an adjacent crystalline cell. 
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I. BACKGROUND 


A. INTRODUCTION 
When a solid is exposed to energetic particles, point 


defects can be produced in the crystal atewee Of the 
material. These point defects can be of various types, such 
as interstitials, substitutional impurities, and vacancies. 
These point defects can interact, and possibly alter the 
properties of the material. This study focuses on a 
molecular dynamics study of vacancies and interstitials in 
the ordered intermetallic compound Zr3Al. This alloy has 
been analyzed in the late 70’s and early 80’s for its 
possible use as a material for pressure tubes in nuclear 
reactors [Ref. 1]. Therefore, studies on the radiation 
damage events and the associated radiation effects on the 
alloy are of particular interest. 

It iS important to point out the difference between 
radiation damage and radiation effects. Radiation damage is 
the term applied to the microscopic events produced by 
energetic particles which irradiate a material. For example, 
a fast neutron may collide with a lattice atom and knock it 
oumeom the lattice, thereby creating a vacancy. The 
displaced atom may take up a position in the lattice that is 
not a regular lattice site. This atom is then termed an 


interstitial. Radiation effects are the macroscopic changes 


of the solid produced as a result of irradiation by energetic 
particles. The radiation damage process is very fast, and 
temperature independent since the event occurs at energies 
much larger than thermal energies [Ref. 2]. Radiation 
effects are temperature dependent, and take a long time to be 
observed. For example, the irradiation of reactor fuel 
elements results in the hardening of the component metals or 
alloys, which leads to an increase in yield strength, but a 
decrease in ductility [Ref. 3]. Under abnormal conditions, 
this property change may be deleterious. However, it is only 


a concern after many hours of reactor operation. 


B. HISTORICAL OVERVIEW 

One of the first studies of the importance of radiation 
effects in solids has its origin in the work begun in 1942 at 
the Metallurgical Laboratories of the Manhattan District in 
Chicago. This laboratory had the task of developing large 
scale reactors which could produce weapons grade plutonium. 
In the middle of 1942, Spedding and Teller [Ref. 4] predicted 
that the high radiation flux generated from the operation of 
large reactors could produce changes in the mechanical 
properties of reactor materials. 

In 1942, Wigner [Ref. 5] did a theoretical study of the 
effects of a large, fast neutron flux on the properties of 
reactor fuel elements. Looking at the fraction of atoms 
displaced from their normal lattice sites in the graphite 


moderator, he concluded that the effects of radiation damage 


Gould not be ignored in“the design and’ operation of large 
scale reactors. From these studies and others, physicists 
and engineers working on the Hanford reactor design were able 
to predict what would happen to the solid components of the 
reactor fuel elements and moderator [Ref. 4]. 

The pioneers in nuclear engineering were not alone in 
studying the effects of radiation damage in solids. The 
irradiation of semi-conductors was also of interest since the 
performance of semi-conductors can be adversely affected 
under bombardment by protons, neutrons, gamma rays and 
electrons. In 1950, Brattain and Pearson [Ref. 6] at Bell 
Laboratories bombarded n-type germanium with alpha particles 
from polonium. A thin slab of material was bombarded on one 
Side and the conductance normal to the bombarding direction 
was measured as a function of time. It was observed that the 
conductance of germanium was Significantly affected under 
irradiation by alpha particles, to the extent that germanium 
converted from n-type to p-type at the minimum conductance 
point. 

Early theories on the mechanism of radiation damage in 
solids were concerned with the production of primary knock-on 
atoms (PKA). In the case of neutron bombardment, the highly 
energetic incident particle undergoes elastic or inelastic 
collisions with the atoms of the "target" material, whereby 
the incident particle transfers enough energy to the target 


aeom townrecmit from adits erystal lattice site. This target 


atom is called the primary knock-on atom, and it can be 
charged, which can lead to the following: 

1. Ionization or excitation en electrons i1nethes@eryseae 

2. The primary knock-on atom can collide with other lattice 
atoms to produce secondary knock-on atoms, which are 
freed from their lattice sites, but with a kinetic 
energy lower than the primary knock-on atom. 

3. The transmission of energy to neighboring lattice atoms 
in an amount which can thermally excite the atom, but 
not enough so as to remove it from its crystal lattice 
Site. 

The third possibility was termed a "thermal spike," since the 
effect of the energy transfer was analogous to 
instantaneously raising the temperature of the atom. In 1949 
Siegel showed the effects of this thermal spike on the order 
in the alloy Cu3Au [Ref. 7]. Siegel showed that the 
irradiation of a solid by the neutron flux in a reactor can 
produce displaced atoms, lattice vacancies, localized thermal 
lattice vibrations, and foreign atoms, all of which can 
affect the macroscopic properties of the alloy. 

These early studies coupled with advances in the fields 
of solid state physics and metallurgy led to more modern 
theories of the mechanism of radiation damage. These 
theories are based upon the defects present in crystal 
Structures, which are formed under irradiation conditions. 

By studying how the defects move and aggregate in the solid, 
a more accurate description of the mechanism of radiation 


damage and subsequent radiation effects can be obtained. For 


example, in 1978, Howe and Rainville [Ref. 8] bombarded 


Z2r3Al, an ordered intermetallic compound, with he Nee and 


Cut ions, and observed the irradiation-induced transition 
sequence ordered — disordered — amorphous state. The 
irradiation-induced crystalline to amorphous transition in 
ordered intermetallic compounds is a widely known phenomenon, 
and much effort is being devoted to understand the mechanism 
Smeets traonsetron (Ref. 99-10]. Thus, it is particularly 
important to understand the point defect configuration in 


these Materials. 


C. PHYSICAL BACKGROUND 

Radiation damage can be defined as the spatial 
distribution of point defects which remains after a solid is 
irradiated. Radiation effects are the macroscopically 
observable changes in the properties and structural features 
of the solid after irradiation of the material. The earliest 
theory of radiation damage was based on a collision theory. 
fies theory, a highly energetic bombarding particle, for 
example a neutron of energy >0.1MeV, collides with the 
nucleus of an atom in the crystal lattice of the target 
material, and the collision results in a transfer of kinetic 
energy to the atom. If the energy transferred is greater 
than the binding energy of the atom in its lattice position, 
the atom will be knocked out of its lattice site and move 
through the crystal. This PKA has an energy on the order of 
Pee eouwlUmereas Of keV [Ref. 11]. The motion of the PKA 
away from its original lattice site results in a vacancy. 


S. 


The PKA will in turn collide with other atoms in the crystal 
lattice and possibly displace other atoms from their original 
lattice sites. Therefore, a cascade o£ atomichecoliaicts es 
the result of a single bombardment event. The displaced 
atoms ultimately appear as interstitial atoms in the lattice. 
A complex type of interstitial is one in the dumbbell 
configuration, which is achieved when two atoms share one 
lattice site. This configuration has been described and 
studied by many researchers [Ref. 12-13]. In summary, a 
Single bombarding event produces an ensemble of point defects 
in the crystal. The subsequent behavior of point defects, 
that is whether they aggregate, annihilate in the material, 
Or are retained, ultimately gives rise to changes in the 
physical properties of the crystal. 

Theoretical approaches to the study of radiation effects 
can involve the concept of a cascade of point defects in the 
crystal. The studies of displacement cascades depend 
critically on the interatomic forces and potentials present. 
In the most simplistic approach, the cascade is the result of 
independent, two-body collisions between the PKA’s and atoms 
in lattice sites. Knock-on atoms are then assumed to move 
freely between collisions [Ref. 14]. This is known as the 
binary collision approximation. A more advanced approach to 
the study of radiation effects is by computer simulations 


[Ref. 15], which take into account not only more realistic 


interatomic forces, but also multiple interactions. This is 
known aS molecular dynamics and is the approach used here. 

It has been found that the type of radiation damage 
observed depends on the material, irradiation temperature, 
and the nature of the irradiating particle. One property of 
materials which is of interest in the nuclear engineering 
field is the amorphization of solids, and under given 
1rradiation conditions, whether or not a material will 
undergo amorphization. 

Amorphization of a crystal is the loss of translational 
symmetry. Some ordered intermetallic compounds which possess 
regular atomic arrangements, a characteristic of a crystal, 
Can undergo amorphization when irradiated. Subsequent 
heating above a certain temperature allows the amorphous 
alloy to crystallize, and radiation damage will have been 
annealed. Annealing is a heat treatment commonly used to 
eliminate structural defects in a metal or to induce phase 
transformations. Table Al gives some examples of ordered 
intermetallic alloys. The ordered intermetallic alloy Zr3zAl 
has been a candidate material for use in nuclear reactors 
because of its low thermal neutron cross section and its good 
tensile properties in the unirradiated state [Ref. 1]. 
Understanding the alloy’s performance in the irradiated state 
tS teus) Of technical importance. 

In 1949, Siegel [Ref. 7] observed that irradiation of 


the ordered alloy Cu3Au can induce the loss of chemical long 


range order if the temperature of the alloy is kept under 


200° CPS but not amorpheyzatioen- The ordered — disordered 
transition was manifested by a reduction in the resistivity 
of an ordered sample after irradiation by fast neutrons to a 
fliuence of 3.3*1019 N/om at 40°CP™ By contrast pen 197 37 
Howe and Rainville [Ref. 16] studied the effects of ion 
bombardment on Z2rzAl. They found that long range order was 
lost. Also, they reported the gradual amorphization of the 
alloy with increasing dose at a rate which depended on the 
temperature of irradiation. 

As previously stated, the point defect configurations 
constitute a step of primary importance in the understanding 
of radiation effects such as amorphization. The reason why 
these configurations are important is because current 
theories of amorphization mechanisms are based upon the 
assumption that a build up of irradiation induced point 
defects is what causes destabilization of the crystalline 


lattice and its collapse to an amorphous state [Ref. 17]. 


D. COMPUTER SIMULATION 

In 1960, Gibson, Goland, Milgram and Vineyard [Ref. 15] 
used a computer simulation to study radiation damage in 
copper. They studied radiation damage events up to energies 
of 400 eV. The code generated the copper lattice and 
Simulated the effects of knock-on atoms produced by 
irradiation. This was the first use of a computer simulation 
for the purpose of studying radiation damage in materials. 
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One of the first attempts at an analytical theory was in 
1954, when Huntington [Ref. 18] developed a theoretical study 
St ehemeacdlatzon damage inJreC metals. One limitation of 
current analytic theory is that the radiation damage process 
is considered as a cascade of hard-sphere, independent 
collisions. Computer Simulations allow one to take into 
account the many-body features of the radiation damage 
process. The high-speed of the computer coupled with its 
large memory allows the characterization of each target atom 
by its mass, interatomic potential, position and velocity. 
This allows Newton’s laws to be solved numerically, as is 
done in the multiple interaction (MI) time-step model used in 
this thesis. The MI code used in this thesis is BCTOOI, 
Peeters a full-lattice simulation developed at the Naval 


POostam~acduate School. 


E. APPLICATIONS 

The study of radiation damage in materials is useful for 
several reasons. It is of fundamental importance with regard 
to the technological industries of nuclear engineering and 
space exploration, but it is also important in the production 
of microelectronic devices. Here materials undergo 
controlled changes after bombardment by ion beams. Ion 
implantation is a controlled technique used in semiconductor 


manufacture. 


Intermetallic phase precipitates in the Zircaloys are 


known to become amorphous after prolonged irradiation 


[Ref. 19]. Zircaloys are very widely used in 
industry. No studies are as yet available on 
these changes on the mechanical properties of 


on the limitations to their use as a function 


the nuclear 
the effects of 
the alloys or 


of dose. The 


Zircaloy matrix undergoes one particular radiation effect 


known aS irradiation growth, a dimensional change which 


limits its useful life in a reactor environment [Ref. 20]. 


The exposure of semi-conductor components 


present in 


Space systems also makes a study of radiation damage in 


materials important. The fluxes of various types of 


radiations which impinge upon the components of satellites or 


manned vehicles can produce harmful effects on semi-conductor 


components. By studying the mechanism of radiation damage in 


semi-conductor components, it may then be possible to develop 


methods, such as the use of shielding, to counter the harmful 


effects of radiation on critical components. 
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Il. OBJECTIVES 


A. THESIS OBJECTIVES 

This thesis is part of a research collaboration program 
initiated by the late Professor D.E. Harrison of the Naval 
Postgraduate School and Dr. D.F. Pedraza of the Oak Ridge 
National Laboratory. The goal of the program is to simulate 
the irradiation-induced crystalline to amorphous transition 
in intermetallic compounds, assuming a mechanism of defect 
buildup. One defect of interest is a vacancy-interstitial 
coupled pair, as described in several papers [Ref. 9-10]. 

In the context of that program, it is first necessary to 
study the relaxed atomic configurations around a vacancy and 
around a cluster of vacancies, as well as the atomic 
displacements involved and the energy of the defects. 
Second, it is necessary to calculate self-interstitial 
configurations, single as well as clusters. Important 
computations include atomic displacements around single and 
clustered self-interstitials, and the corresponding energies 
Prcdmmoremeteon volumes... Third» it is necessary to calculate 
the vacancy-interstitial spontaneous recombination distance 
as a function of the self-interstitial orientation. The 
Orientation of concern is the dumbbell configuration, with 
the distance mentioned above measured relative to the 
Ve@eamey. This study will deal with the first and third 


Prepeems . 


ame 


B. PREVIOUS EFFORTS 

Although there has been extensive work reported on the 
Simulation of point defect configurations and the 
displacement field around them, most has been done using a 
Static lattice approach [Ref. 21]. In this approach, the 
defect is "moved" into different positions until the 
potential energy of the lattice reaches a minimum. Although 
there have also been molecular dynamics simulations of point 
defects, most studies use "static" approaches, where a 
minimum of the lattice potential energy is sought, 
irrespective of the dynamics by which that minimum can be 
achieved. In these cases, the Simulation progresses until 
the total force acting on every atom of the crystallite is 
zero. A few studies involving the use of such methods are 
described below. 

In 1968, Englert, Tompa and Bullough [Ref. 22] studied 
point defects and dislocations in copper. They developed a 
pair potential for copper which was used to study the atomic 
configurations of various line and point defects in the 
metal. A crystallite containing the defect of interest was 
Simulated, for FCC copper. The crystallite was allowed to 
find its minimum potential energy configuration via a 
relaxation procedure known as the conjugate gradient 
technique [Ref. 23]. This technique involves the evaluation 
of the first derivative of the potential, the force, and is 


designed to guarantee the convergence of the potential energy 


2 


to a minimum energy. They found that the minimum energy 
atomic configuration under conditions of a single vacancy 
corresponded to 1.09eV. This was in good agreement with the 
experimental value of 1.14eV. They also computed the atomic 
displacement field around the vacancy, up to the seventh 
nearest neighbor, which is contained in Table A2. 

Tn 1974, Savino and Perrin [Ref. 24] studied the 
morphology of planar vacancy aggregates in copper. The 
computer simulation created a FCC lattice with the vacancy 
defect of interest, and the method of conjugate gradients was 
used to obtain the equilibrium energy configuration. This 
study found that six vacancies collapsed to form either loops 
or stacking fault tetrahedra. 

The method of conjugate gradients has also been used to 
Study the effects of interstitials on the atomic 
Comenaquratton and energies in FCC metals. In 1976 Schober 
[Ref. 13] studied single and multiple interstitials in 
copper, using the computer program DEVIL developed in AERE 
Harwell. Computed values included the atomic displacements 
and relaxation volumes for various defect configurations. 
Various results are shown in Table A3. Single and multiple 
interstitials were found to be highly mobile, in that the 


activation energies for migration of the defect are low. 


C. A DYNAMIC SIMULATION 
The studies above are "static" approaches, since a 


minimum potential energy of the lattice was sought. The 


3 


simulation is Garried out untiieehe total torce Germngeon 
every atom in the crystal is zero. In this thesis, a fully 
dynamic simulation is attempted. The defect is introduced 
into the crystallite, resulting ina perturbed crystal. As 
the simulation progresses, the perturbation decays and the 
atoms are displaced to positions in which they are fairly 
stable. 

For the relaxed vacancy studies, the perturbation 
involves removing one or several atoms from their normal 
lattice sites. The atomic displacements of surrounding atoms 
are then computed. For the studies of vacancy-interstitial 
recombination, the perturbation involves placing an 
interstitial in the dumbbell configuration in a crystal 
lattice, with a corresponding vacancy. One atom of the 
dumbbell is given a kinetic energy, and this perturbation and 
its effects on the lattice are be studied. Of concern is 
the spontaneous recombination distance as a function of the 
dumbbell orientation relative to the vacancy. 

Most simulations of point defects have been done for pure 
metals, or metals containing small amounts of impurities. 

The reason is that for fully quantitative studies, an 
accurate interatomic potential is necessary. In this study, 
the interest is in intermetallic compounds, which are 
concentrated alloys. Since this study will not approach the 


more difficult subject of obtaining accurate potentials for 


14 


alloys, the focus will be on qualitative and semi- 


quantitative results. 
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Pil. COMPUTER SIMULATION AND MODEL DEVELOPMENT 


A. THE COMPUTER MODEL AND RELATED PROGRAMS 
Lew eCT 

The computer simulation used in this study is called 
BCT, which is a program developed by the late Professor D. E. 
Harrison of the Naval Postgraduate School. This program used 
multiple-interaction (MI) timestep logic, since many events 
occur simultaneously. The program BCT allows the 
Simultaneous computation of the position and velocity of an 
atomic particle, which are obtained via a numerical solution 
to Newton’s laws of classical mechanics [Ref. 25]. 

The initial input data to the program vary slightly, 
depending on the topic of study. For vacancy relaxation 
studies, the initial input data include the target 
crystalline structure, atomic masses, elements to be included 
in the crystal, and the locations of the vacancies. For the 
spontaneous-recombination studies, additional input includes 
the initial energy provided to one atom of the dumbbell 
interstitial, as well as the relative position of the 
interstitial in the crystal lattice. The atom which receives 
the momentum impulse may travel in any direction in the 
lattice. The momentum is input using two angles, theta and 
phi, which are one direction angles in a spherical coordinate 


System. In the code, theta is the angle between the momentum 


IS 


vector and the positive y-direction, and phi is the angle 
Beeweetmeene pPpOsittive x-Girection and the projection of the 
momentum on the surface plane. The geometry is indicated in 
Peguce Bl. 

The computer program also has the capability of 
including inelastic energy loss. In all studies, the input 
variable ICVIN will determine whether or not inelastic losses 
due to electronic effects are to be included. The Lindhard 
Inelastic Loss Model [Ref. 26] is used to simulate electronic 


energy losses, which arise due to interactions between moving 


charges and free electrons in the lattice. The losses are 
represented in terms of the electronic stopping power, S,, 
which is defined as: 
See wena ae za zeee. (1) 
where an icucohemBboha radius, 
ape nh“ /ameme* |; (2) 
and, 
e = the electronic charge 
m= the electron mass 
Z, = atomic number of a novuengeparticle 
Z5 = atomic number of the lattice atom 
h = Planck’s constant 
ZS ait eile me peal 3/2 ey 


The constant Zo¢¢ is the effective charge of the two 
interacting atoms. For all studies, the value for S, was 


1.44*10723 kg/s. The electronic stopping power is 


1a 


essentially a friction force, and acts to dampen the forces 
that arise due to the interatomic potentials. The model 
assumes that the electronic stopping power is proportional to 
the velocity of the moving atom [Ref. 27]. In all further 
discussions, electronic energy losses will be referred to as 
EE losses. 

The program develops the collision cascade while 
continuously tracking the positions and velocities of lattice 
atoms in time. If the positions and velocities were computed 
for all lattice atoms every timestep, excessive computer run 
time would be used. Therefore, forces are not computed on a 
lattice atom unless this atom is within 2.4 LU from a moving 
atom. The program maintains a list of neighbors to moving 
atoms, which is periodically updated. In addition, to save 
computer time, the timestep increment is variable, and is 
obtained by dividing a specified distance by the speed of the 
highest atomic velocity in the crystallite. The specified 
distance is 0.1 LU, with the lattice unit being defined as 
one-half the lattice parameter, ag: The program 1s run for a 
specified number of timesteps, usually enough to establish a 
lattice in equilibrium, or to infer the approach to 
equilibrium. 

2. Ancillary Programs 

Several other subprograms are used in the simulation 

to support the main program. Several of the important 


programs are discussed below. 
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a. INTEG 

This program performs all of the numerical 
integrations encoded in the simulation. Newton’s laws of 
motion are used to calculate the positions and velocities of 
atoms in the lattice. 

beeePOTVEEC 

This program calculates the forces and potential 
energies for all moving atoms and their neighbors. Forces 
are obtained by interpolating force tables, and potential 
energies by integrating force tables. 

es POTIN and POTFEN (RK) 

These programs involve the identification of 
input potential function parameters, and the creation of the 
potential tables. The forces on atoms in the crystallite are 
obtained from these tables. 

digest BETOO 

This program generates a body centered tetragonal 
(BCT) lattice, which is entirely equivalent to an FCC 
structure. The crystallite faces are (100) pee = (010) Rap, 
(Oll) pee = (001) pep, and (011) peq= (100) per: 

e. HALLOY 

This program generates an Ll, crystal structure 
as indicated in Figure B2, when an ordered lattice is 
selected. When a disordered lattice is desired, Al atoms are 
inserted at random locations in the crystal, keeping the 


composition Zr7cAloc. 


i) 


B. PROPERTIES OF Zr AND Zr3Al 

A summary of various properties and physical constants 
for Zr and Zr3Al is shown in Table A4 which was obtained from 
several sources [Ref. 28-30]. 

Table A3 shows that there are two equilibrium crystalline 
Structures for pure Zr, hexagonal close packed (HCP) below 
860°C, and body centered cubic (BCC) above this temperature. 
The phase diagram for the Zr-Al binary system is shown in 
Figure B3. The beta phase corresponds to the HCP structure 
of Zr, and the alpha phase corresponds to the BCC structure. 
The figure also shows that Z2r3Al forms via a peritectoid 
reaction [Ref. 28], 

Z2r5Al + Zr => Z2r3zAl 
A peritectoid reaction is one in which two solids combine to 
produce a third solid on cooling. The peritectoid 
equilibrium is shown to occur at about 9/5 Gaewiened 
corresponding value of about 9% Al by weight. 

This study involves the use of the program BCT001 which 
generates a body centered tetragonal crystal equivalent to 
the FCC structure which corresponds to ZrzAl. This is 
evident in Figure B4, which contains the trigonometric 


relationships between the FCC lattice and the BCT lattice. 


C. POTENTIAL FUNCTIONS 
1. Basics 
Multiple Interaction computer simulations involve the 
computations of forces between atoms. The forces are 
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obtained from a potential table, which in turn is formed by 
evaluating a chosen potential function. Therefore, if 
realistic data are desired, one must have an accurate 
potential function for the system being studied. 
Unfortunately, the exact form of interatomic potential 
functions and forces between two atoms are not known. The 
reason 1s that the problem is a complex, many-bodied problen, 
in that there are interactions between electrons with their 
respective nuclei, between electrons of one atom, and between 
electrons of both atoms. The problem is resolved by making 
approximations and assumptions, which lead to the selection 
of one of the many theoretical and empirical interatomic 
potentials available. Many different types of potential 
functions have been utilized in radiation damage studies. 
The potential functions used in this study are discussed 
briefly below. More detailed descriptions of the potential 
functions can be found in Torrens [Ref. 31]. 

a. Moliere 

The Moliere potential function is a purely 

repulsive potential function. It is convenient to use this 
function to represent the wall of an interatomic potential. 
It is categorized as a "screened coulomb potential," and is 
an approximation to the Thomas-Fermi screening function 
(Ref. 31]. Basically, it is used for small interatomic 


distances, whereby the nuclei of the two atoms exert 


aM 


coulombic forces on each other. The form of the screened 
interatomic potential is: 
Vie (Gri <a aii (4) 

The general form of the Moliere interatomic potential is: 

V(r) = (Z,25e%/r) (0.35exp (-0.3xr/a) + (5) 

0.55exp (-1.2r/a) + O.lexp (-6.0r/a). 

The variable "a" is the Firsov screening length, which is 
defined as: 

a = k ON85Sap/Me! =e eee (6) 
and usually k=1, but here k will be an adjustable variable 
near 1 for providing a smooth spline, to the potential well. 

b. MORSE 
The Morse potential function is both an 
attractive and repulsive function, depending on the 
interatomic distance. This function is used in sputtering 
Simulations and radiation damage studies as the attractive 
part of the interatomic potential function. It has the form: 
V(r) = D exp[-2a(r-R,) ] - 2D, exp i agai CSR aa (7) 


e 
D, 1s the depth of thespereneratwely Ra Sei equilibrium 
interatomic separation distance, and ais usually the bulk 
modulus or the compressibility. However, Harrison [Ref. 32] 
finds it more beneficial to leave "a" as an adjustable 
parameter. The parameter "a" can then be used to match this 


attractive Morse potential to the repulSive Moliere potential 


at some separation. 
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c. Composite Morse-Moliere 
In this study, the repulsive wall of the Moliere 
function is joined with one cubic spline to the attractive 
Moers] function. First, the a term in the Morse and the k 
term of the Moliere are adjusted to allow a smooth match 
between the two functions, and then the cubic spline connects 
the two. In this way, there exists a potential function over 
a finite range of interatomic separation distances, and 
molecular dynamics studies of vacancies and interstitials can 
me Conducted . 
2. Selection of Potential Function Parameters 

A composite Morse-Moliere interatomic potential 
function was used in this study. The potential function used 
Pere simulation trials is indicated in Figure BS. The 
potential function used was the same, regardless of whether 
the code was run using pure Zr, or with the alloy. This will 
be discussed further below. The potential parameters chosen 
are indicated in Table A5. The terms a and k were selected 
i onmeer CO Obtain a smooth potential function curve. The 
well depth value was chosen to match the cohesive energy for 
purer er, which 1S 6.316 eV. The equilibrium separation 
distance was set to 3.09 A (the lattice parameter of Zr3Al1 is 
aWee2en). The values R, and R,, are the spline boundaries. 
Between these values the spline matches the two curves. The 
value R, 1s a cut-off distance beyond which potential and 


POrees are Zero. 
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IV. RESULTS 


A. RELAXED ATOMIC CONFIGURATIONS 
1. General 

Various vacancy sites were used in order to study the 
relaxed atomic configurations around a vacancy or vacancy 
cluster. This study focuses on single vacancies, 
di-vacancies, and tri-vacancies. Various input parameters 
were changed in order to study their effects on the 
displacements of atoms surrounding the vacancy or cluster of 
vacancies. For example, simulations were performed with and 
without frictiongtorces, sand witheanwalle, Someries a... 

In each trial performed, the displacements of atoms 
Surrounding the vacancy or vacancies were computed. Average 
atomic displacements were computed for first nearest 
neighbors and second nearest neighbors. Nearest neighbors 
are lattice atoms separated by 1.414 LU, and second nearest 
neighbors are separated by 2 LU. In each tri-vacancy run, 
the nearest neighbors of one vacancy were identified so as to 
track the distance of each neighbor atom to the vacancy each 
timestep. Plots were made of this separation distance versus 
the elapsed time. This was done in order to visualize the 
effects of various parameters on the oscillatory motion of 


the lattice atoms. 
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2. Single Vacancy Results 

One configuration was used to study the atomic 
displacements caused by a single vacancy. The atom which was 
removed from the lattice to create a vacancy was the body 
Semeenea atom Of the BCT unit cell. This configuration is 
portrayed in Figure B6. Two trials were performed, as 
indicated in Table A6é. Trials 1 and 2 were performed with no 
EE losses. Trial 1 was performed with a Zr lattice, and 
trial 2 was performed with an alloy. 

The quantity of interest in these two trials was the 
effect of the crystal type on the average atomic 


displacements, which are produced when the lattice relaxes to 


amiemestablemeontveuratwonm. Specifically, the difference in 
mass, or the "mass effect," between the two runs, was of 
interest. Since the interaction potential between the Zr and 


Al atoms are the same, the static configuration of the alloy 
should be the same as that of the pure metal. However, in 
inducing the motion of atoms by the creation of a vacancy, 
then the differing mass of the two particles will affect 
their dynamics. This could possibly have the effect of 
producing different configurations for the relaxed crystal, 
both of which correspond to local minima in the potential 
energy. 

The spatial arrangement of atoms about the single 
vacancy iS portrayed in Figure B6. The body centered 


position has 12 nearest neighbors, and 6 second nearest 
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neighbors. As shown in Figure Bo, there are two classes of 
first nearest neighbors, which are marked Fl or F2. The 
class Fl corresponds to first nearest neighbors which occupy 
body centered positions in cells adjacent to the one 
containing the vacancy. There are 4 class Fl first nearest 
neighbors. The class F2 corresponds to the 8 atoms which lie 
on the corners of the BCT unit cell containing the vacancy. 
There are 6 second nearest neighbors comprising class Sl. 
The reason for this classification scheme is because it is 
found that the neighbors in these sites have different 
relaxed configurations whose final displacements are 
approximately the same. 

The average atomic displacements for class Fl, F2, 
and Sl are contained in Table A6. For Zr, or Trial 1, the 
first nearest neighbors comprising class Fl show the largest 
average atomic displacement, -0.0159 LU, which is inward, or 
towards the vacancy, aS indicated by the minus sign. The 
first nearest neighbors comprising class F2 show an inward 
displacement of -0.0074 LU. Thus, a first nearest neighbor 
occupying a body centered position, adjacent to a cell 
containing a vacancy in a body centered position, has the 
larger displacement. The second nearest neighbors, class Sl, 
show an outward average displacement of +0.0065 LU. These 
results are in qualitative agreement with the results of 
copper [Ref. 22]. For trial 2, performed with an alloy, the 


average atomic displacements are, in general, smaller than 
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for Zr. However in both cases the atoms in the relaxed 
configurations had lower potential energies than before 
relaxation. For example, in both the alloy and Zr lattice, 
the value for the initial potential energy of S2 atoms was 
-5.73 eV. Typical values for the potential energy after 
PeMeratwonwwere —=5981 tow=5e91 eV for the alloy, and from 
—weeetnor —5.71 eVefor Zr. 
3. Di-Vacancy Results 

Two distinct configurations were used to study the 

atomic displacements caused by a di-vacancy. The vacancies 


were located at sites marked with open symbols in Figures B7 


and B8. Figures B7 and B8 correspond to configurations Dl 
and D2, respectively. Six trials were performed, three on 
each configuration. In two of the trials, the run was 


performed with EE losses, but with either an alloy or pure 
Zr. This was done to determine if there was a noticeable 
effect of the mass difference under conditions of EE losses. 
PaleQ,e for eacheconfigumation, two trials with pure Zr, one 
with and one without EE losses, were compared. This was done 
to see if there was an effect of EE losses on the final 
atomic displacements. 

The initial input parameters and average atomic 
displacements for di-vacancy configuration Dl are contained 
in Table A7. In this configuration, there were five classes 
of nearest neighbors, Fl, F2, F3, Sl and S2. Examples of the 


spatial arrangement of four of these classes are contained in 
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Figure B7. The atoms labelled Fl in Figure B7 correspond to 
first nearest neighbors of one vacancy. These atoms lie in 
a body centered position of a unit cell adjacent to the 
vacancy, as shown in Figure B7. There are eight class Fl 
atoms. The atoms labelled F2 are first nearest neighbors to 
one vacancy, which lie in the corners of adjacent unit cells. 
There are eight class F2 atoms. Atoms labelled F3 are first 
nearest neighbors to both vacancies. These four atoms lie in 
the body centered plane located between the two vacancies. 
Atoms labelled Sl are second nearest neighbors, and lie in 
corners of adjacent unit cells. There are eight class Sl 
atoms. Finally, atoms listed as S2 in Table A7 are second 
nearest neighbors to each vacancy, which are found by moving 
2 LU in the +y om -y direction from the pos@telonmer eacn 
vacancy. There are two class S2 atoms. 

Table A7 shows that for each trial run, the average 
atomic displacement obtained for each nearest neighbor class 
did not significantly differ. The average atomic 
displacements were dependent on the spatial location of the 
atom about the vacancy. For example, Fl atoms were the only 
class which showed an inward displacement, that is, towards 
the vacancy. The value obtained for the average atomic 
displacement was on the order of -0.016 LU. Class F2 atoms 
moved away from the vacancy, on the order of +0.006 LU. 
Class F3 atoms moved away from the vacancy, on the order of 


+0.002 LU. Class Si atoms moved away from the vacancy, on 
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Pos Ordenr@ot Peo wOolL 29 LUM Panally, class S2 atoms"moved away 
Epemethe vacancy, on the order of +0.004 LU. 

The input parameters and average atomic displacements 
for di-vacancy configuration D2 are contained in Table A8. 

In this configuration, there were five classes of nearest 
Wemegnbors, F1,8F2, 9S, Sl and S2. Examples of the spatial 
arrangement of these classes are contained in Figure B8. The 
atoms labelled Fl in Figure B8 correspond to the body 
centered, first nearest neighbors of one vacancy. There are 
twelve class Fl atoms. The atoms labelled F2 are first 
nearest neighbors to one vacancy, which lie in corners of 
adjacent unit cells. There are four class F2 atoms. Atoms 
labelled F3 are also first nearest neighbors to one vacancy, 
an@ware located by moving 1 LU in the +z or -z direction from 
each vacancy. There are two class F3 atoms. Atoms labelled 
Sl correspond to second nearest neighbors, and lie in corners 
of adjacent unit cells. There are four class Sl atoms. 
meal, atvomsslabelled S2 lie an corners of unit cells, 
Wiehwacesfound by moving 2 LU in thesty or -y direction from 
Pacmmacancy. There are four class,S2 atoms. 

Table A8 shows that for each trial run, the average 
atomic displacements obtained for each nearest neighbor class 
did not significantly differ. This was the same result 
obtained from configuration Dl trials. Again, the average 
atomic displacements were dependent on the spatial location 


of the atom about the vacancy. For configuration D2, the 
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body centered, first nearest neighbors comprising class Fl 
moved inward, on the order of -0.017 LU. All other neighbors 
moved outward, and class F2 atoms moved away from the vacancy 
as much as class Fl atoms moved towards it. 

In general, for di-vacancy studies, atoms located in 
body centered positions of cells adjacent to the cell 
containing the vacancy show the largest displacements. In 
both configurations, the direction of the displacements was 
towards the vacancy. The results also did not show any 
Significant effect of EE losses or mass on the average atomic 
displacements for each class of nearest neighbors. 

4. Tri-Vacancy Results 

Two distinct configurations were used to study the 
atomic displacements caused by a tri-vacancy. The vacancies 
were located at sites marked with open symbols and a label 
"vy" in Figures B9 and B10. Figures B9 and B10 correspond to 
configurations Tl and T2, respectively. Eight trials were 
performed, four on each configuration. For each 
configuration, two trials were performed with Zr, one trial 
with EE losses and the other without. The other two trials 
were performed with the alloy, one case with EE losses, and 
the other without. In this way, for each tri-vacancy 
configuration, investigations were made for observing any 
effects of EE losses or mass on the average atomic 


displacements. 
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The input parameters and average atomic displacements 
iu iri—vacane yw conriguration Tl are contained in Table AQ. 
In this configuration, there were four classes of nearest 
neighbors that were considered, Fl, F2, F3 and Sl. Examples 
of the spatial arrangement of these classes are contained in 
Figure B9. The atoms labelled Fl in Figure B9 occupy body 
centered positions of the Six unit cells adjacent to the cell 
Gemearning™the tri-vacancy. The six atoms*of class Fl are 
all first nearest neighbors to either vacancy V, or V2. The 
atoms labelled F2 are first nearest neighbors to one of the 
corner vacancies. There are Six class F2 atoms. Atoms 
labelled F3 are the eight nearest neighbors to the body 
centered vacancy alone. Finally, atoms labelled Sl are 
second nearest neighbors to one of the vacancies. There are 
eight class Sl atoms. 

Table A9 shows that for most runs, the average atomic 
displacements obtained for each nearest neighbor class did 
not significantly differ. The average atomic displacement 
was dependent on the spatial location of the atom about the 
vacancy. The only trial which did not follow this trend was 
trial 3. In trials 1, 2 and 3, the largest average atomic 
displacements were directed inward, towards the vacancy. 
First nearest neighbors of class F2 and second nearest 
neighbors of class Sl show outward displacements, with a 


magnitude less than that for class Fl. 
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The input parameters and average atomic displacements 
for tri-vacancy configuration T2 are contained in Table Al10O. 
In this configuration, three classes of nearest neighbors 
were considered, Fl, F2 and Sl. Examples of the spatial 
arrangement of these classes are contained in Figure B10. 

The atoms labelled Fl in Figure B10 occupy body centered 
positions of cells adjacent to vacancy V1 and vacancy V3. 
There are eight class Fl atoms. The atoms labelled F2 occupy 
body centered positions of cells adjacent to vacancy V2. 
There are four class F2 atoms. Thus, class Fl and class F2 


atoms are first nearest neighbors to the corresponding 


vacancies. Finally, atoms labelled Sl are second nearest 
neighbors to either of the vacancies. There are 14 class Sl 
atoms. 


Table A100 shows that the average atomic displacements 
obtained for each nearest neighbor class did not 
Significantly differ. The one exception was trial 3, where 
the average atomic displacements were low. The reason for 
this discrepancy has not been determined. The average atomic 
displacement was dependent on the spatial location of the 
atom about the vacancy. Atoms of class Fl moved towards the 
vacancy, on the order of -0.012 LU. Atoms of class F2 also 
moved towards the vacancy, by about -0.005 LU. Second 
nearest neighbors of class Sl moved away from their 


respective vacancies, by about +0.0115 LU. 
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5. Atomic Oscillations 

In order to illustrate the transition to 
equilibrium, plots such as those indicated in Figure Bll were 
obtained. This was done for each tri-vacancy run. A 
reference point, the body centered vacancy of the tri-vacancy 
defect, was selected. Several first nearest neighbors 
associated with this body centered vacancy were chosen, and 
the position of the atom from the reference point was plotted 
versus time. In this way, equilibrium could be seen by a 
dampening of the oscillations of the atom about its 
equilibrium position. This was readily apparent in trials 
performed with EE loss. For non-EE loss trials, the 
oscillations continued for a longer time until equilibrium 
was reached, as 1S expected. As an example of these plots, 
Figure B12 shows the plot of a class F2 first nearest 
neighbor, and its distance from the body centered vacancy of 
configuration T2. The data was obtained from a lattice of 
pure Zr and no EE losses. Figure Bl2 shows that equilibrium 
occurred at about 10 picoseconds. The large reduction in 
oscillation amplitude define the equilibrium point. Figure 
B13 shows the same atom tracked with EE losses incorporated. 
In this case, equilibrium occurred in about 3 picoseconds. 

In addition to determining the point of equilibrium, 
plots as described above were generated to investigate if 
EWere was an erfect of mass on the oscillatory motion of 


atoms about the vacancy. MTri-vacancy configuration T2 
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trials were compared. One trial was performed using an 
alloy, and one with Zr. For both cases, no EE losses were 
used. Three atoms, which corresponded to geometrically 
distinct lattice sites about the vacancy, were selected for 
observation. For the alloy, the corresponding Figures are 
B14, B15 and B16. For Zr, the corresponding Figures are B17, 
B18 and B19. The periods of the atomic oscillations about 
the equilibrium position are contained in Table All. The 
periods do not Significantly differ, and are between 0.17 and 
0.2 picoseconds. These are of the same order of magnitude as 
typical phonon frequencies although none have been explicitly 


determined for Zr. 


B. RECOMBINATION STUDIES 
1. Lattice Orientations and Terminology 

Two lattice orientations were used in order to study 
vacancy-interstitial spontaneous recombination. The first 
orientation is indicated in Figure B20, and will be referred 
to aS case 1. The second orientation is shown in Figure B21, 
and will be referred to as case 2. For both cases, the 
dumbbell consisting of the interstitial pair is centered 
about the body centered position of the BCT unit cell, with 
its axis parallel to the y-direction. The separation 
distance between the two atoms of the dumbbell was 0.7LU. 
The difference between the two caseS was in the placement of 


the vacancy. 
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FOr@@ase |] @erials, the vacancy was located in a 
corner of the BCT unit cell, which is a nearest neighbor 
position to the body centered position. For case 2 trials, 
the vacancy was located in the body centered position of a 
cell adjacent to that of the dumbbell, as indicated in 
Figure B21. In order to be consistent in discussing the 
atoms of the dumbbell, the one with the largest y-coordinate 
is referred to as the lead atom, which is marked with an "L" 
iperigures B20 and B21. 

In order to visualize the results, two types of 
graphics were used. Three dimensional plots containing the 
initial and final positions of the vacancy’s neighbor atoms 
were constructed. An example is shown in Figure B22. The 
imveial positions are solad circles, and the final positions 
are marked with arrowheads. Vacancies are indicated with 
open circles or squares. When necessary for discussion, Al 
atoms for alloy trials are appropriately indicated. The 
other graphics output used was a plot of the separation 
distance versus elapsed time. The separation distance is 
defined as the distance between the lead atom of the dumbbell 
Aicene Init val position of the vacancy. 

2. Trial Description 

In short, the lead atom of the dumbbell is given an 
energy impulse of a specific magnitude. The values of theta 
and phi allow the lead atom to follow any desired direction. 


Betas yare performed for a lattice consisting of the 2Zr-Al 
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alloy or pure Zr, with and without friction siorcec mas 
vacancy is positioned in one of two positions, corresponding 
to case 1 and case 2 trials. 

A summary of the trials performed along with the 
initial conditions is listed in Table Al2. The column headed 
momentum direction was obtained from selecting various theta 
and phi values. For example, a value of 0 for both theta and 
phi implies the momentum was in the +ty-direction. 

3. Recombination Results 

Table Al3 contains a list of the trials performed as 
well as the elapsed time to equilibrium, and the approximate 
separation distance at equilibrium. In recombination 
studies, equilibrium is defined as the time for the 
separation distance to be 0.15 LU or less, with small 
oscillations about this position. The value chosen is 
arbitrary, and serves as a basis for comparisons between 
various trials. 

Trials performed under conditions in which only one 
input variable differed were compared. The results are 
described below. 

a. Momentum Direction 

In order to determine the effect of the initial 
momentum direction on recombination, comparisons were made on 
trials with the same initial configuration. The only 
variable in the runs was the momentum direction. For 


example, case 1 EE loss runs performed with an energy impulse 
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Gem. l CVgwemomcompaarea. These correspondseto trials 2,4,6 
am@ezion In all cases, the time required for equilibrium was 
2.3 picoseconds, and the separation distance was 0.11 LU. 
mroure B23 shows the in@tial and final positions of the 
dumbbell and neighbors for trial 2, and comparisons of the 
peemepenms Lometwials 4,6 ande26 didmot sagnificantly differ. 

Similar comparisons were made for case 1 trials 
10,12 and 14. In this case, the time required for 
equilibrium was 2.3 picoseconds, and the separation distance 
between the vacancy and interstitial was 0.11 LU. Figure B24 
shows the initial and final positions of the dumbbell and 
Weagnibors for trial 14, and comparisons of the positions for 
metals LO and’ ha) Gad net seaigqnaficantiky differ. 

Thus, it appears that for the given trial 
configuration, the final position of the lead atom is 
independent of its initial direction, and hence its inertia. 
The simulations were terminated after 1200 timesteps, usually 
around 5 picoseconds. Recombination occurred to the extent 
that there was always some small separation between one atom 
of the dumbbell and the vacancy. However, relaxation towards 
equilibrium might either take many more timesteps or never 
Coceus. This Ms aneintrinsac aspect of the molecular dynamics 
approach since a large set of atoms have been set in motion 
and need to disperse the initial imported kinetic energy. 
"Perfect" recombination between the vacancy and dumbbell 


cannot be seen during the timescale of the simulation. In 
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all trials performed, the second atom of the dumbbell 
repositioned itself to the body centered position of the BCT 
unit celle 

De eOrlentation 

In order to determine the effect of defect 
orientation on recombination, comparisons were made between 
case 1 and case 2 trials. 

Trials 1 and 20 were non-EE loss trials. 

Figure B25 shows the initial and final positions of trial 1 
at time 2.94 picoseconds, and Figure B26 shows the same for 
trial 20 at time 3.27 picoseconds. Both "snapshots" of the 
lattice were taken at approximately equal times. The most 
noticeable difference between the two is that the 
recombination is in amore advanced state for trial 20. From 
Figure B25, recombination is still occurring. The separation 
distance between the lead atom in the dumbbell and the 
vacancy is 0.25 LU. Figure B25 also shows much motion for 
neighboring atoms. Thus, there seems to be a preferred 
Orientation for recombination, that being the dumbbell as 
oriented as in case 2. 

Trials 2 and 21 were trials conducted with EE 
loss incorporated in the simulation. Figure B21 shows the 
initial and final positions of trial 2 at time 3.92 
picoseconds, and Figure B27 shows the lattice for trial 21, 
at time 4.34 picoseconds. Again, recombination occurs 


quicker for trial 21. For trial 2), equa briunveceursese 
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time 0.80 picoseconds, after which there are small 
oscillations of the lead atom about the position of the 
Meaeeawey. BPOr the case 1 trial, equilibrium occurred at 2.3 
picoseconds, with a separation distance of 0.11 LU. 

Trials 5 and 23 were non-EE loss trials. 

Figure B28 shows the initial and final positions of trial 5 
at time 3.47 picoseconds, and Figure B29 shows the same for 
trial 23, at time 2.92 picoseconds. Again, it can be seen 
that the preferred orientation for recombination is a case 2 
pesmeton. “Equilibrium for trial 23,.occurred at 0.9 
picoseconds, with a separation distance of 0.07 LU. The 
separation distance for trial 5 was 0.21 LU. 

©: Welectronic Energy Losses 

In order to determine the effect of EE losses on 
recombination, comparisons were made between trials using the 
same lattice orientation, energy, and momentum direction. 
Several pairs of data sets were analyzed. The pairs below 
were selected because they were done with different 
Orientations and energies. This was done to show that the 
general effects of EE loss on recombination were the same 
regardless of other input parameters. 

Trials 17 and 18 were case 1 trials with the lead 
atom of the dumbbell possessing an energy impulse of 0.001 eV 
Peeetem-arrect10On. Trial 17 is the non-EE loss case, and 
trial 18 the EE loss run. Recombination is more likely to 


occur in the EB loss run than in the non-EE loss run. For 


SY, 


the EE loss case, the equilibrium separation distance was 
0.11 LU, which occurred at time 2.3 picoseconds. In the non- 
EE loss case, the equilibrium separation distance was 

0.22 LU. Figures B30 and B31 show the lattice positions for 
the non-EE loss and EE loss case, respectively. The lead 
atom of the dumbbell travels closer to the vacancy site in 
the friction case than the non-EE loss case. Also, the other 
dumbbell atom repositions to the body centered positions more 
evidently in the EE loss case. These effects probably arise 
due to the increased range of motion achieved in the absence 
of EE loss forces. Atoms will travel further and reach 
closer interatomic distances, and hence higher forces, in 
nNONn—-EEe less Gemiars. 

Trials 20 and 21 were case 2 trials with an 
energy impulse of 0.1 eV in the y-direction. In both cases, 
recombination occurred at 0.9 picoseconds, but in the EE loss 
case, the lead atom of the dumbbell is closer to the vacancy 
Site by 0.04 LU in that time. 

In general, trials performed with EE losses 
result in a more defined recombination within the time of the 
Simulation. EE losses also reduce the oscillations of atoms 
about their equilibrium sites. 

di ee oy, 

In order to determine the effect of the mass of a 

lattice atom on recombination, various trials were performed 


using an alloy. Comparisons were made on trials which 
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differed only in the selection of either an alloy or pure Zr 
lattice. Several trials were chosen for discussion. 

Trials 16 and 36 were case 1 trials with the lead 
atom of the dumbbell possessing an energy of 0.01 eV in the 
y-direction. For the alloy run, trial 16, the time to 
equilibrium was 0.5 picoseconds. For pure Zr, equilibrium 
@eeurred at 1.0 picoseconds. The separation distance for the 
alloy and pure Zr was 0.04 LU in each case. Figures B32 and 
Peeeere illustrations of the atomic positions for the case of 
the alloy and Zr lattice, respectively. As seen from these 
plots, it does not appear the use of the alloy had an effect 
emeene recombination process. 

Trials 21 and 33 were case 2 trials with the lead 
atom of the dumbbell possessing an energy of 0.1 eV in the y- 
Gaeectuen. AS indicated in Table A100, equilibrium in both 
cases occurred at the same time and at the same separation 
distance. Figures B28 and B34 show the lattice positions, 
[heGae support this observation. 

Mregeneral, 2 seems that the use of an alloy or 
pure Zr does not significantly affect the recombination of 


mae Vacancy-interstitial defect. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The results of this study show that the molecular 
dynamics approach is a valid tool which can be utilized in 
the study of atomic relaxation and vacancy-interstitial 
recombination. 

Comparisons of the atomic relaxation studies show that 
the final average atomic displacements are not significantly 
affected by inelastic energy losses or the introduction of a 
mass effect. Rather, the atomic displacements are highly 
dependent on the initial orientation of the atom with respect 
to the vacancy. In general, atoms that are first nearest 
neighbors to a vacancy, and which occupy the body centered 
positions of cells adjacent to the vacancy, will move the 
most, and will move inward, towards the vacancy. Second 
nearest neighbors generally move away from the vacancy. 
These results are in qualitative agreement with static 
relaxation studies performed with single vacancies in FCC 
copper. 

The studies performed on vacancy-interstitial 
recombination show that recombination is not affected by the 
inertia or the mass effect of the alloy. Electronic energy 
losses result in a more defined recombination during the 
timescale of the simulation. Also, there does seem to be a 


preferred orientation of the vacancy-interstitial defect, for 
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the two cases investigated. Recombination occurs faster with 
the vacancy oriented in a body centered position in a cell 
adjacent to the dumbbell. This is supported by the findings 
of the relaxation studies. The atoms which occupy the body 
centered positions adjacent to a body centered vacancy move 
the most. Since the dumbbell is oriented about the body 
Gemteerecasp@e@sittion Of a Unit cell, and has a large potential 
energy due to the small interatomic distance, 0.7 LU, it is 
expected that the atoms of the dumbbell will move quickly to 
minimize the energy of the system. 

Given that the atomic displacements obtained in vacancy 
relaxation studies are symmetrical, and correspond 
G@@alitatively with others results, this molecular dynamics 
technique should be used to study defect configurations. The 
code could be used for computing the average atomic 
displacements for lattices containing multiple vacancies. 
Also, the distance between the vacancy and the dumbbell could 
be manipulated in order to establish the distance beyond 


which recombination does not occur. 
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APPENDIX A - TABLES 


TABLE Al. ORDERED INTERMETALLIC COMPOUNDS 


Intermetallic Compounds Gry Sta lSseruct ume 


NiAl B2 
Na nal Ll 
Co la B2 
Co5T1 Calls 
Cra ia CALS 
ema Bil 
FeTi B2 
Neo a: B2 
Nia EST E93 
Nacsa DO 
NWA cee 
2r5Al ll 
Z2r5Al Ll, 
Ue e D2, 


for copper is 3.603 A. 
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TABLE AZ DISPLACEMENTS OF ATOMS IN SHELLS AROUND A SINGLE 
VACANCY IN FCC COPPER. 
Shell 
First S.6 @On ce Manis FOUGER. .-.Lithe ote Seventh 
-0.019 +0.002 -0.003 -0.005 0 -0.002 -0.001 
Note: Units are in terms of the lattice constant, which 


TABLE A3. 


DISPLACEMENTS DUE TO A <100> DUMBBELL WITH A 


MODIFIED MORSE POTENTIAL 


shell 


dumbbell 


O11 
110 


Note: Displacements 


for copper 


TABLE A4. 


Atomic Number 


Atomic Weight 


Density 
Cryoua! 


Lattice 


Valence 


(g/cm?) 


Type 


Constant 


(2) 


(on 


(amu) 


(A) 


Displacement (a) 
S Ss = 


<<< —_  _ yy 
0 0 O'S 
0 O08 ORO 

-0.03 -0.03 8) 


are in terms of the lattice constant 
603 A). 


PHYSICAL DATA FOR Zr AND Z2r3Al 


Zr 2r3Al 
40 
Cale 22 
ie 2o5 Be 0 
rer, BCC FCC (ordered) 
aj=3.231 (HCP) Aceon 2 


Coy -144 (HCP ) 
a_j=3 . 62 (BCC) 


Hy) tue 
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TABLE AS. POTENTIAL PARAMETERS 


Equilibrium 
Well Separation Adjustable Spline Cue—oOf £ 
Depth Distance Parameters Boundaries Distance 
(eV) (A) ae) (LU) (LU) 
D, Re a k a Ry Sos 
Ge eoc 308 gS 0.9 Ol fshs 0.92 2.4 


TABLE A6. AVERAGE ATOMIC DISPLACEMENTS FOR SINGLE VACANCY 
RELAXATION STUDIES 


Average Atomic Displacements for 
Various Nearest Neighbor Classes 


(LU) 
EE 
Trial Crystal Losses F1 (4) F2 (8) S1 (6) 
ak Wise no -.0159 -.0074 00.6.5 
2 Alloy no -.0004 ta O10 tT 00Ze 
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TABLE A7. AVERAGE ATOMIC DISPLACEMENTS FOR DI-VACANCY 
RELAXATION STUDIES (CONFIGURATION D1) 


Average Atomic Displacements for 
Various Nearest Neighbor Classes 


(LU) 
EE 
ieral Crystal Losses F1(8) F2 (8) F3 (4) S1 (8) SB) 
it Zr no Oeste) “teO027 +40124 +.004 
Z Zid yes Pe oeeteOdg? “4.0020 +.0118 +.004 
3 Alloy yes = OLS) @, OOo?) =. O023° +20120 +.004 


Note: The number in parentheses next to the nearest neighbor 
class represents the number of atoms in the class. 


TABLE A8. AVERAGE ATOMIC DISPLACEMENTS FOR DI-VACANCY 
RELAXATION STUDIES (CONFIGURATION D2) 


Average Atomic Displacements for 
Various Nearest Neighbor Classes 


(LU) 
EE 
Meial Crystal Losses F1(12) F2 (4) Poe ) S1 (4) S2 (4) 
Ak Ze no =e? Pteebos +2,0061 +.0718 =+.004 
2 ar yes =sol73s teureo +.0064 +.0074 +.004 
3 Alloy yes ta set aooe +0065 +,0094 +.004 


Note: The number in parentheses next to the nearest neighbor 
class represents the number of atoms in the class. 
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TABLE A9. AVERAGE ATOMIC DISPLACEMENTS FOR TRI-VACANCY 
RELAXATION STUDIES (CONFIGURATION T1) 


Average Atomic Displacements for 
Various Nearest Neighbor Classes 


Vi) 
EE 
Trial Crystal Losses F1(6) F2 (6) F3 (8) S1 (8) 
it ax no —.0252 +.0168 ~+2007 0s, - CZ 
2 Jae yes -.0152 +.0131 +.003 +.0264 
5 Alloy eke) +.0046 +.0195 -.0140 -.0045 
4 Alloy yes =-.0152 +.0131 +.0030 +.0133 


Note: The number in parentheses next to the nearest neighbor 
class represents the number of atoms in the class. 


TABLE A10Q. AVERAGE ATOMIC DISPLACEMENTS FOR TRI-VACANCY 


RELAXATION STUDIES (CONFIGURATION T2) 


Average Atomic Displacements for 
Various Nearest Neighbor Classes 


(LU) 
EE 
Trial Crystal Losses F1(8) F2 (4) SsiL(r4) 
il Ze no =e Olea =. O11 S se Ol Ik yi 
2 Ze yes = 0b 23 = ONOONS + Ales 
& Alloy no -.0066 -.0012 +.0108 
4 Alloy yes en) eZ SOOO IG ee 2 


Note: The number in parentheses next to the nearest neighbor 
class represents the number of atoms in the class. 
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PAB All. Oscillations 


Period (picoseconds) 


Atom Alloy Ar 
i Crd eo, 18 
2 220 uals 
S O23 @= 20 
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TABLE Ald2. 
SlUDIES 


Trial 


WMATA BWNHE 


Gryseal 


Alloy 
Alloy 
Alloy 
Alloy 
Alloy 


ALISON | 
ALLOy 


Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 


AEOy® 
NLL 


Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
Alloy 
ar 
2X 
Ze 
ZX 
ie 
4x 
Zs 
Ae 


indicates a disordered alloy 


Friction 


no 
yes 
no 
yes 
no 
yes 
no 
yes 
no 
yes 
no 
yes 
no 
yes 
no 
yes 
no 
yes 
no 
no 
yes 
no 
no 
no 
no 
yes 
yes 
yes 
yes 
no 
yes 
no 
yes 
yes 
yes 
yes 


Case 


PROMO RPRFRPRPMONMENNNNHNNFPRPRP ERP REP EP RP RP PRP RP PPP PPP PS 
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Energy 
(eV) 


Ose Sete Oo 2 2 2 2 OO Oro OO ©. 2 2 0 OO OO © OC Cree] 2 CoO 


Theta 


180 


180 


0 


INITIAL INPUT PARAMETERS FOR RECOMBINATION 


Momentum Direction 


Phi 


0 ME eo A De E> | a A a JE as al a ED 


TABLE A1l3. EQUILIBRIUM TIME AND SEPARATION DISTANCE FOR 
RECOMBINATION STUDIES 


Time to Reach Equilibrium Separation 
Equador eum Distance 
irial (picoseconds) (LU) 
a N = 
2 ag O.241) 
3 N = 
4 OE real I 
5 N = 
6 23 Spanien 
7 3 CralZ 
8 O25 Olas 
9 
10 285 Orel 
1 N = 
£2 ZS Ora 
i 3 N = 
14 255 0 apie 
a5 0.5 020s 
16 ORS 0.04 
i] N = 
18 265 Ogi 
19 N = 
20 ieee) 02.06 
21 0.8 oe02 
AL 0.6 0.06 
23 O79 b707 
24 Ores Crus 
25 O38 02.02 
26 Pe O71 
an Oe: 0.03 
28 digs | 0.04 
29 One 0.04 
30 eg 6. 10 
ou 170 0.04 
o2 dees oz 10 
33 G9 C01 
34 ae C202 
Se 1.0 0.04 
36 ee) 0.04 


Seb 
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The Direction Angles Theta and Phi in Spherical 
The Plane is the Surface Normal 


Figure Bl. 
The Direction y is into 
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Figure B2. Illustration of a L12 (Ordered FCC) Lattice. The 
Length a, is the Lattice Constant. 
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The Zirconium Potential Function. 
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Pigure@Bs6:. 


The Configuration for Single Vacancy Atomic 


Relaxation Studies. The Vacancy is Marked v 
and Examples of Atoms of Nearest Neighbor Class 


Fl, F2 and 81 are Illustrated. 
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Figure. B7. 
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The Spatial Arrangement of Atoms for Di-Vacancy 
Configuration Dl. The Two Vacancies are Marked 
V and Examples of Atoms of Nearest Neighbor 
Class Fl, F2, 


F3 and 81 are Illustrated. 
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Figure B8. 


The Spatial Arrangement of Atoms for Di-Vacancy 
Configuration D2. The Two Vacancies are Marked 
V and Examples of Atoms of Nearest Neighbor 
Class Fil, F2, F3, 81 and 82 are Illustrated. 


5 2) 





Figure B9. The Spatial Arrangement of Atoms for Tri-Vacancy 
Configuration Tl. The Three Vacancies are Marked 
V1, V2 and V3. Examples of Atoms of Nearest 
Neighbor Class Fl, F2, F3 and 81 are Illustrated. 
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Figure Bio. 


, (3 





The Spatial Arrangement of Atoms for Tri-Vacancy 
Configuration T2. The Three Vacancies are 
Marked V1, V2 and V3. 


Examples of Atoms of 
Nearest Neighbor Class Fl, F2 and 81 are 
Illustrated. 
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Figure Bll. Sample Plot of the Separation Distance Between 
the Neighbor Atom and Body Centered Vacancy. 
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Elapsed Time (picoseconds) 


The Oscillation of a Class F2 Nearest Neighbor 
About the Body Centered Vacancy of Configuration 
T2. The Plot Was Obtained from a Lattice of 
Pure Zr Without EE Losses. 
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Figure Bl3. The Oscillation of a Class F2 Nearest Neighbor 
About the Body Centered Vacancy of Configuration 
T2. The Plot Was Obtained from a Lattice of 
Pure Zr With EE Losses. 
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Figure B15. 
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Oscillation of Atom 2 About the Body Centered 
Vacancy. The Plot Was Obtained Using the Alloy 
With No EE Losses. 
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Figure B16. Oscillation of Atom 3 About the Body Centered 
Vacancy. The Plot Was Obtained Using the Alloy 
With No EE Losses. 
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Oscillation of Atom 1 About the Body Centered 
Vacancy. The Plot Was Obtained Using Zr With 
No EE Losses. 
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Figure B18. Oscillation of Atom 2 About the Body Centered 
Vacancy. The Plot Was Obtained Using Zr With 
No EE Losses. 
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Oscillation of Atom 3 About the Body Centered 
Vacancy. The Plot Was Obtained Using Zr With 
No EE Losses. 
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Figure B20. Vacancy-Interstitial Configuration for 
Recombination Case 1 Trials. 
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Vacancy-Interstitial Configuration for 
Recombination Case 2 Trials. 
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Figure B22. 


Sample Plot of the Atomic Displacements 


Produced as a Result of Vacancy-Interstitial 
Recombination. 
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Figure B23. Atomic Displacements for Recombination Trial 2. 
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Figure B24. Atomic Displacements for Recombination Trial 14. 
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Figure B25. Atomic Displacements for Recombination Trial 1. 


76 


7 POSITION (LU) 





Figure B26. Atomic Displacements for Recombination Trial 20. 
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Figure B27. Atomic Displacements for Recombination Trial 21. 
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Figure B28. Atomic Displacements for Recombination Trial 5. 
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Figure B29. Atomic Displacements for Recombination Trial 23. 
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Figure B30. Atomic Displacements for Recombination Trial 17. 
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Pigure B32. 





Atomic Displacements for Recombination Trial 16. 
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Figure B33. Atomic Displacements for Recombination Trial 36. 
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Figure B34. Atomic Displacements for Recombination Trial 33. 
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